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SUMMARY 


A 0.137^n NACA 0012-64 airfoil has been tested in the Langley 0.3-m 
transonic cryogenic tunnel at free-stream Mach numbers of 0.75» O. 85 , and 0.95 
over a total pressure range from 1.2 to 5.0 atmospheres. The onset of conden- 
sation effects were found to correlate more with the amount of supercooling 

in the free-stream than it did with the supercooling in the region of maximum 

\ 

local Mach number over the airfoil. Effects in the pressure distribution over 
the airfoil Were generally seen to appear over its entire length at nearly the 
same total temperature. Both observations suggest the possibility of hetero- 
geneous nucleation occurring in the free-stream. Comparisons of the present 
onset results are made to calculations by Sivier and data by Goglia. The 
potential operational benefits of the supercooling realized are presented 
in terms of increased Reynolds number capability at a given tunnel total 
pressure, reduced drive fan power if Reynolds number is held constant, and 
reduced liquid nitrogen consumption if Reynolds number is again constant. 
Depending on total pressure and free-stream Mach number, these three benefits 
are found to respectively vary from 7 to 19 %, 11 to 25 %, and 9 to 20 $. 
Appendices are included which give details of the data analysis procedure 
and of the error bar calculation for the differences in pressure distributions. 



INTRODUCTION 


Cryogenic wind tunnels should he operated at the lowest possible total 
temperatures in order to maximize the advantages of cryogenic operation, such 
as increased Reynolds number and reduced operating costs. However, minimum 
operating temperatures are limited by low temperature effects on the flow, , 
which can be the result of real-gas behavior of the test-gas. Either the equa- 
tion of state of the test-gas changes in such a manner that the gas does not 
properly simulate the nearly ideal-gas behavior of air encountered in flight, 
or the gas begins to condense and consequently does- not properly simulate the 
nearly ideal-gas behavior of air. Analyzing. the nitrogen test-gas over the 
operating ranges of existing and proposed transonic cryogenic wind tunnels, 
Adcock reported in reference 1 that low temperature nitrogen gas approximates 
an ideal-gas during flow simulations for total pressures under 10 atmospheres 
and for Mach numbers up to 2. Consequently, the lower temperature boundary 
appears to be determined by condensation of the nitrogen test-gas. The 
boundary is important from an operational viewpoint because it determines, on 
one hand, the maximum Reynolds number capability for a fixed tunnel total 
pressure and, on the other hand, the most economical total temperature and 
pressure to operate the tunnel for a fixed test Reynolds number. 

Condensation not only has an important role in determining the operational 
envelope for transonic cryogenic tunnels but has a similar role for supersonic - 
and hypersonic tunnels. Wegener and Mack in reference 2 give a summary of 
historical events leading to, and the results from, many of the condensation 
investigations in supersonic and hypersonic tunnels. 
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The present condensation investigation employs a pressure-instrumented 
NACA 0012-61; airfoil mounted in the Langley 0.3-m transonic cryogenic tunnel. 
Changes in the pressure distribution over the airfoil have been used to detect 
the onset of condensation effects in a similar fashion as changes in pressure 
distributions in divergent nozzles were used in references 2 to 6 to detect 
the onset of condensation effects at supersonic and hypersonic speeds. The 
airfoil has been tested at free-stream Mach numbers of 0.75, 0.85, and 0.95 
over a total pressure range of 1.2 to 5*0 atmospheres in order to determine the 
total temperature at which condensation influences the flow about the airfoil 
as a function of total pressure and free-stream Mach number. 
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sound, speed 

airfoil chord, 0.137-m 

P ~ Poo 

pressure coefficient , 

difference in C defined “by equation (A6) 
P 

pressure, atm 

pressure instrumentation accuracy 
dynamic pressure 
Reynolds number 
temperature , K 

supercooling defined by equation (l) 
linear dimension along airfoil chord line 
viscosity 
density 

standard deviation 


Subscripts 

e conditions at onset of condensation effects 

L local conditions 

s . conditions at intersection of isentrope and vapor pressure curve 

t total conditions 

00 free— stream conditions 



MINIMUM OPERATING TEMPERATURES 


The maximum Reynolds number capability and the operating costs of a 

cryogenic wind tunnel are both affected by the minimum operating temperatures 

of the tunnel. The reasons for this can be understood by reviewing the effects 

of changing tunnel temperature upon the test gas, which are shown in fig- . 

ures 1(a) and l(b). In figure l(a) the behavior of the gas properties relative 

to their ambient temperature values can be seen. The density, p, increases 

while the viscosity, y, and the speed of sound, a, decrease as the total - 

temperature, T^., decreases for a fixed value of tunnel total pressure, p^. 

As shown in figure 1(b), the behavior of the gas properties as T^ decreases 

leads to a dramatic increase in unit Reynolds number, R, a decrease in drive 

power required, and no change in dynamic pressure, q. Consequently, it is 

desirable to operate the tunnel as cold as is possible in order to maximize R 

and reduce drive power. Just how low in temperature the tunnel can be operated 

is determined by the onset of condensation effects. 

A conservative .first guess for the minimum operating temperature is to 

operate at temperatures which avoid any possibility of saturation occurring 

« 

anywhere over the model. Since the lowest static temperature over a model 

occurs at the point of maximum local Mach number, M^ , T and p can be 

max 

chosen to keep the local static values of p, T just on the vapor phase side 
of the vapor pressure curve for nitrogen. However, it has been found from tests 
in the Langley 0.3-m transonic cryogenic tunnel and in other facilities that 
condensation effects tend to be delayed to some lower temperature. 

The delay in condensation onset past the vapor— pressure curve is due to 
an energy barrier blocking stable liquid droplet formation - that is, blocking 
self-nucleation of the test gas. Thus there will normally be a delay in 
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droplet formation and growth, and consequently in effects, unless some catalyst 
is available to speed the process. (More details of the various condensation 
processes will be given in the section Implications of the Data. ) The delay 
in condensation can be quantitatively measured by the supercooling at onset of 
effects, AT, which is defined as 

AT = T - T (1) 

s e 

'where T is the static temperature vhere the isentrope crosses the vapor 
s 

pressure curve and T g is the static temperature along the isentrope where 
effects are first seen. 

If supercooling does occur then there are two wind tunnel operational 
benefits that result. First of all, for a given. tunnel p^ the maximum R 

capability is .increased because at the lower temperature the unit . R is greater 

as shown in figure l(b). This benefit would not only increase the maximum R 
capability of a tunnel, but would also increase the R obtainable for a model 
that may be q limited since q is proportional to p. for a given M . 
Secondly, if maximum R capability is not important, then supercooling can be 
used to reduce operating costs when operating at a fixed R. Operating at a 
lower Tj. than needed to avoid saturation permits a decrease in p^ for a 

given value of R. Since drive fan power is proportional to p^ and decreases 

with decreasing T^_ , both changes reduce the required drive power. Since 
normally over 99 % of the injected liquid nitrogen is used to absorb the heat 
generated by the drive fan once the tunnel has been cooled, a reduction in 
drive fan power leads to a reduction in the amount of injected liquid nitrogen. 
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Both reductions decrease direct operating costs, which are split with about 
90$ going to the injected liquid nitrogen and only 10$ going to power the 
drive fan. Specific details of drive power required and operating costs are 
given in references T and 8. 
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EXPERIMENTAL APPARATUS 


Tunnel 

The Langley 0.3-m transonic cryogenic tunnel is a continuous flow, fan- 
driven tunnel, which uses nitrogen as a test gas and is cooled by injecting 
liquid nitrogen directly into the stream. Variation of the rate of liquid . 
nitrogen injection provides a total temperature range from nearly 77 to 350 K, 
while the total pressure can he varied from 1.2 to 5*0 atm. The combined low 
temperature and high pressure can produce a Reynolds number of over 330 million/m 
(.100 million/ft). Some of the design features and operational characteristics 
of the 0.3-m tunnel have been reported by Kilgore in reference 9, and a sketch 
of the tunnel is shown in figure 2. 

The liquid nitrogen used to cool down the tunnel and absorb the heat of 

\ 

the drive fan was injected into the tunnel through a series of eleven nozzles 
arranged on three different struts, or spray bars, at the three injection 
stations shown in figure 2. Full details of this arrangement are included in 
reference 9* However, for one of the tests during the present investigation 
(details to be mentioned in RESULTS AND DISCUSSION), the spray bar system was 
removed and injection was carried out by injecting through just four nozzles 
at injection station 1. This change, of course, necessitated larger liquid 
flow rates per nozzle. 


Airfoil and Installation 

A 0.137-m NACA 0012-64 airfoil was used for these tests. Starting at the 
leading edge, the airfoil had 20 pressure orifices spaced at 3 % chord intervals 
on both the top and bottom surfaces. For the tests at of 0.75 and 0.95, 

a rearward-facing orifice was added to the trailing edge of the airfoil. The 
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airfoil was installed "between flats in the octagonal test section, with the 
leading edge 0.62-m from, the "beginning of the test section. The angle of 
attack of the airfoil was zero for all tests. A sketch of the airfoil is 
given in figure 3. 

The pressure distributions over the airfoil may be influenced by wall 
interference effects because of the relatively large chord of the airfoil 
compared to the size of the test section. However, wall interference should 
not have an effect on the occurrence of condensation. 

Data Acquisition and Error Discussion 
The tunnel parameters were first recorded and then a differential pressure 
transducer was used in a scanning valve system to step through and record the 
pressures around the airfoil. The time to acquire all of the data for a 
complete pressure distribution was 50 seconds. 

The uncertainty in the pressure transducer measurements was less than 0.5 percent 

of full scale, which was the magnitude specified by the manufacturer. In fact, the seatte 

in the pressure transducer data was observed to be less than the stated value by more 

than a factor of 10. For the transducer used in the tests at M m = 0.75 and 

0.95» the uncertainty was on the order of 0.0004 atm. For the transducer used in 

the test at = 0.85, the uncertainty was about 0.0002 atm. There was no 

significant error introduced by either the signal conditioning or data acquisition 

systems. However, during the 50-second acquisition period, the tunnel conditions 

were observed to fluctuate by the following amounts: M^, ±.003; T^, ±0.5K; and 

P+s approximately ±0.005 atm. A development of the uncertainty in C due to the 
u p 

above fluctuations is contained in appendix B. 
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TESTS 


Data Sampled 

To determine the total temperature at which effects do occur and to inves- 
tigate the influence of total pressure and free-stream Mach number on the onset 
of condensation effects, the airfoil was tested at- free-stream Mach, numbers of 
0.75, 0.85, and 0.95 over a total pressure range of 1.2 to 5.0 atm. For each 
of the three values of test envelopes’ were drawn which covered those 

tunnel p^ and, ranges that were of interest and that were obtainable in 
the 0 . 3 -m tunnel. 

The test envelopes were bounded in p^ by the 1.2 to 5.0 atm pressure 
limits of the tunnel and were bounded in T^ by total temperatures corre- 
sponding to saturated flow at and those corresponding to saturation 

max 

in the reservoir section. Figure 4 shows the M - O .85 test envelope as an 

example. Total conditions which fall lower in than the reservoir 

saturation line 'are unobtainable since the injected liquid nitrogen used for cooling 

would never evaporate when T is below the vapor pressure curve and, consequently, furtt 

cooling can not occur. Total temperatures above those along the line of local 
saturation offer no possibility of condensation in a pure nitrogen test-gas so 
they are used in this study only as a way of obtaining unaffected comparison 
points. The locus of tunnel p^ and T^ which just allow saturation when 
the flow is accelerated to is also given in figure 4 . 

To explore each of the regions of interest in the and plane, 

paths of constant R and M were used to traverse the region of interest. 

Since along each of the paths both M ro and R q were held constant," systematic 
deviation in pressure coefficients at any x/c position over the airfoil was 
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.taken to -"be the result of condensation effects. The paths of constant R 

c 

,and .are shown for the = 0.85 test in figure 5. The R c along each 
path yaried from l6 million at the low pressure traverse to k2 million at the 
high pressure traverse. Along these 5 paths, .a total of 201 pressure .distri- 
butions -were recorded and analyzed. Similar envelopes were drawn ;and 
investigated for .= 0.75 and = 0.95 except that only 3 paths ,of 
.constant *R and M were used in these cases. These envelopes are shown in 

, - » C CO ' , - 

figures ,6 and J. During the investigation of the high pressure path for the 
M m = 0.95 test, a fan drive power limit was encountered so M ot was reduced 
from 0.95 to 0.93 for that path alone. 

Data Analysis 

In order to determine the p^ and at which the onset of condensation 

effects occurs along a constant and R c path, a special data analysis 

procedure was developed. The main objective of the data analysis procedure was 

to be able to correct values of C on an orifice-by-orifice basis in cases 

P 

where the pressure distribution over the airfoil was taken at a slightly 
different than desired. With such a correction, data along the path that 

normally would have had to be disregarded due to small differences in M could 

CO 

be utilized. Having more data available made it easier to detect the total 

temperature at which differences in C , AC , begin to occur because of the 

P P 

effects of condensation. 

Once, the values of total temperature for the onset of effects at each 
orifice had been determined, a graph of orifice onset total temperature as a 
function of x/c. could be drawn for the path of' constant M and R . In 

03 c 

this. manner,, an appropriate onset total temperature could be chosen for the 

t 

airfoil as- a whole. Further details of the data analysis procedure are 
presented in appendix A. 
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RESULTS AND DISCUSSION 
Data Presentation and Comments 

The data, taken and analyzed as discussed in the previous section, -will be 
summarized by a series of figures for each value of - 0.75, 0.'85, and 0.95- 
covered in this study. Each series will include average, condensation-free 
pressure distributions for each path of constant R c ; graphs of orifice onset 
Tj. as a function of x/c for each path; and a summary plot for each 
showing the onset conditions for each path in a T versus p plane. The 
results for all of the M ro will then be grouped together assuming either 

max 

or M o 3 -*- s important Mach number parameter. 

The first test described is the M^ = 0.75 test. The average pressure 

distributions are shown in figures 8, 9s and 10, and the orifice onset T, vs. 

t 

x/c graphs are shown in figures 11, 12, and 13. Four different symbols are 
plotted in the orifice onset graphs corresponding to four different situations. 
When normal orifice onset is detected, a circle is plotted at the appropriate 
x/c, Tj_ location. When no orifice, onset is detected, a square is plotted at 
the minimum T^ at which data was taken. When excessive scatter obscures any 
trends, a triangle is plotted at the minimum T, . When an orifice pressure 

■ U 

tube developed a leak, a diamond is plotted at the minimum T . It is 

"G 

interesting to note in figures 11, 12, and 13 that the airfoil appears to be 
most sensitive to condensation effects in the 15-25# chord location, as is 
seen most clearly in figures 12 and 13. There does not seem to be any greater 
sensitivity to condensation over the rearward portion of the airfoil, which 
seems to downplay the importance of any condensate growth over the airfoil 
itself. The summary plot of onset conditions ds shown in figure l4. 


13 



The pressure distributions for the 5 lines of constant R c for the 
M ot = O.85 test are shown in figures 15 , l6, IT, 18, and 19 and exhibit a 
recompression shock at the O.65 x/c location. The graphs showing orifice 
onset T^ as a function of x/c are shown in figures 20, 21, 22, 23, and 24. 
With the possible exception of figure 24, it appears that effects 
occur at about the same T. over the length of the airfoil. That is, there is 
not ah obvious sensitivity to condensation effects behind the shock, or even 
right before the shock in the region of high local Mach number. Along the. 

R c = 42 million path as shown in figure 24, there does appear to be possible 
sensitivity behind the shock location. At x/c values of 0.75 and O.85 the 
onset value of T^ appeared to be 100.5 K, which is above the trend of effects 
occurring at about 96 or 96.5 K for the majority of orifices. While the onset 

T^ of 100.5 K may be correct, there is scatter in the AC^ vs. T^ graphs that 

were combined to form figure 24. Because of the uncertainties involved and 
because of the lack of agreement between orifice onset Tj. ' s behind the shock 
location at x/c =0.65 in figure 24, the airfoil onset T was chosen to be 
98.5 K with a possible error bar extending from 96.5 to 100.5 K. The summary 
plot of onset conditions for the = 0.85 test is given in figure 25. 

The last test condition was that of M m = 0.95. As mentioned under the 
section Data Sampled, the high pressure path at R c = 40 million was tested 

at M^ = 0.93 instead of 0.95, so there is a difference between the first 

2 paths and the R c = 40 million path. The unaffected pressure distributions 
are shown in figures 26, 27, and 28, and the orifice onset T^ vs, x/c graphs 
are presented' in figures 29, 30, and 31. ' I!he onset T t 's for the R q = 16 million 
path, do not seem to have any one area more sensitive than any others, as 




14 



•seen in. figure 29 . However, onset T’s for the R - 2T million path in 

t c 

figure 30 may be showing signs of an increased sensitivity at about x/c — O.iO 
in the region of increasing local Mach number, although the trend 

is not conclusive. The third path at R c = 40 million shows a very uniform 

onset T^ over the airfoil. There was no data taken along this path between 

T = 100.5 and T = 103.5; consequently, there is no resolution of effects 

between these 2 temperatures and this may contribute to the fairly uniform 

distribution of onset T^. The lack of resolution is indeed unfortunate as 

seen in figure 32, which shows AC vs. T, for x/c = 0.95, because there 

P x 

are no effects at T^_ - 103.5 K but well developed effects at 100. 5 K. For 

this path, onset T^. for the airfoil was chosen to be 102.0 K with the error 

bar spanning 100.5 to 103.5 K. The summary plot for onset conditions for the 

M» = 0.95 test is given in figure 33. 

The results for all of the tests will now be grouped together assuming 

first that is the important Mach number parameter ana then assuming 

max 

that is the important parameter. The question of which Mach number is 

important arises because of the uncertainty in the mode of nucleation that 'is 
occurring, which is discussed in the next section. 

The first pair of figures comparing all of the onset results is constructed 
for the supercooling as defined in equation (l). 

AT = T - T (1) 

s e ■ 

In equation (l), the value of T g is the static temperature at which the 

isentrope, starting at the onset values of p^ and T^., crosses the vapor- 

pressure curve. The value of T g is the static temperature along the same 

isentrope at M = in figure 34(a) and at M = in figure 34(b). 

max 
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In both cases' the value of AT is plotted as a function of the onset P along 

each path of constant R c . In figure 34(a), the supercooling was nearly constant 

at 7 K over the pressure range for M m = 0.75, was roughly equal to 11 K over 

the pressure range for = 0.85, and was 12 and 14 K for the two lower 

pressure paths at = 0.95 and then dropped to 10 K for the high pressure path. 

In figure 3Mb), the supercooling based on M^ remained between 2 and 3 K over 

the pressure range for both = 0.75 and' O.85. However, the supercooling for 

the M^ = 0.95 test was 1 and 2 K for the lower pressure paths, but dropped to 

—2 K for the upper pressure path. Regardless of the behavior of the high pressure 

path for the = 0.95 test, there seems to be much better correlation of the 

values for supercooling when assuming M^ is the pertinent Mach number rather 

than M^ 

max 

The second pair of figures again compares all of the onset results, but the 

comparison is now made in a log^Q p versus T plane. Static values of p, T 

are calculated from the p^_ and T^ at onset and are calculated assuming an 

isentropic expansion to either or M^. Figure 35(a) shows the p, T 

max 

as calculated assuming M, is the’ Mach number of interest while figure 35(b) 

max 

shows the p, T as calculated assuming M^ is the important Mach number. It 

would again appear from figure 35, as in figure 34, that the onset data correlates 

much better with M^ than with M^ 

max 

The seemingly early onset of effects in the high pressure path of the 
M^ = 0.95 test may be the result of a tunnel configuration change that was made 
before the data for this path was taken. The spray bars system used for injection 
of the liquid nitrogen as described under the tunnel apparatus was removed for 
another investigation in the .3-m tunnel and remained out for the tunnel entry 
during which the data was obtained for the high pressure path of the M^ = 0.95 
test.. With the spray bars removed the liquid was injected from four nozzles 


16 



flush, with the interior tunnel wall at injection station 1 instead of the 
eleven nozzles at the three different stations. Consequently, the liquid 
injected may not have "broken up and evaporated as effectively as it did with 
the eleven nozzles. The unevaporated drops may have formed the increased number 
of nucleating sites necessary to decrease the amount of supercooling realized. 
Because of this uncertainty due to the tunnel configuration change, the data 
from the high pressure path of the = 0.95 test will -not he presented in 
further data comparison graphs. 

Implications of the Data 

Condensation effects can he caused by three different modes of nucleation, 
where nucleation is taken as the formation of a stable liquid droplet. The 
three different modes are homogeneous, heterogeneous, and binary nucleation. 
These three processes differ by the manner in which they overcome an energy 
barrier to droplet formation and hence growth. In depth discussion of the 
three types can be found in references 2, 3, 10, 11, 12, and 13. The 

different modes can briefly be summarized as follows. 

Homogeneous nucleation can occur when the gas undergoes supercooling, which 
will decrease the energy barrier to liquid droplet formation. When enough 
supercooling occurs, the gas molecules will be able to overcome the energy 
harrier to droplet formation and massive nucleation rates are likely. Hetero- 
geneous nucleation may occur when there are liquid or solid impurities 
suspended in the vapor at the time the flow first expands past the vapor- 
pressure curve and begins to supercool. The energy barrier of formation, 
which is so important for homogeneous condensation, is not so significant here 
because the impurities can act as a catalyst for the condensation process. 
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Consequently, little or no supercooling is needed before stable droplet growth, 
occurs. What prevents heterogeneous nucleation from always dominating the 
homogeneous process is that it requires extremely large concentrations of 
impurity sites to influence the flow as was pointed out by Wegener and Mack 
in reference 2 . The third type of nucleation possible is binary. This process 
occurs when a mixture of two different gases condenses into droplets which are 
composed of both types of molecules. The binary process is discussed 
specifically in references 10 to 13. The combination of the two species can 
effectively change the energy barrier to nucleation in such a manner that 
condensation will occur earlier than analysis would predict for either of the 
species acting singly. 

The data from the present experiment does seem to suggest which nucleation 

process might have occurred. First, no condensation seems to be occurring until 

at least 7 K of supercooling has taken place at the location. If binary 

max 

nucleation is occurring with some sort of impurity vapor in the nitrogen used 
for cooling the tunnel, then its effects were either undetectable op occurred 
at low enough temperature to be confused with homogeneous or heterogeneous 
nucleation of the nitrogen test gas. The purity of the liquid nitrogen injected 
for cooling was approximately 99-95 percent in the present experiments, so it is 
assumed that binary nucleation would be unlikely. Second, none of the orifice 
onset T^ versus x/c plots show any systematic trend toward condensation first 
appearing in the region of maximum supercooling, which is the region of high 
local Mach number. Since the amount of supercooling is so important for 
homogeneous nucleation, these results do not seem consistent with what one would ' 

t 

expect if homogeneous nucleation was occurring. In fact, there seems to be no 
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generally preferred location over the airfoil for effects to occur. Hence the 
effects may actually be the result of heterogeneous nucleation occurring upstream 
of the model in the test section. Under these circumstances the airfoil would 
experience the influence of these droplets over its length at about the same 
although local sensitivity to the condensate in the free— stream may vary over the 
airfoil because of differences in or pressure gradient. If condensate in 

the free-stream is the cause of the effects, the possible source for the impurity 
sites could be the injected liquid nitrogen used for cooling the tunnel. At the 
very low temperatures of these tests, the injected liquid may not be completely 
evaporating in the time it takes to travel from the injection stations to the 
test section (see figure 2). The partially evaporated droplets may provide the 
high concentration of nucleation sites needed for heterogeneous condensation. 

Some evidence for heterogeneous nucleation from unevaporated droplets in the 
free-stream comes from figures 3^ and 35 » which show stronger correlation : 
between the onset of effects and than the onset of effects and 

max 

Further evidence comes from the smaller amount of supercooling realized with 
the high pressure path of the = 0.95 path. Because the configuration without 
spray bars is less efficient at breaking up the injected liquid nitrogen, it 
would seem that the smaller amount of supercooling was due to an increased 
number of liquid drops on which condensation growth could occur. Light 
scattering tests should be made in the test— section to determine conclusively 
if indeed there are unevaporated injection droplets. 

Comparison to Other Works 

While it appears that heterogeneous nucleation is the cause of condensation 
in the present experiments, most of the previous condensation investigations have 



19 



concentrated on homogeneous nucleation. Some of these 'investigations were 

concerned with nitrogen gas, and two references, l4 and 15, will "be used here 

for comparison to the present data. A traditional means of presenting onset of 

effects data is in the log.^ P vs. T plane. As mentioned in the previous 

section, there is some question as to the process of nucleation and hence in 

the relevant condition to plot in the log_^ Q p vs. T plane. For the present 

comparisons, it was decided to use the static values of p and T as calculated 

for- the M^ condition "because regardless of the source of the effects,- this 
max 

amount of supercooling was realized in the region of maximum local Mach number 
before effects were detected. 

i 

Using the log^ Q P vs. T format, the present results can be compared with 
the analytical nitrogen-gas work done by Sivier in reference l4. His computer 
computations assumed homogeneous nucleation and used classical liquid droplet 
theory.. An unpublished curve fit by Adcock to Sivier' s calculated onset points 
for conditions of higher pressure and ‘lower onset Mach- number has been plotted 
along with the present, onset conditions in figure 36, where Adcock’s curve- 'fit 
is represented by 

7jT* = 1.^61 x 10 2 - 4.,g62 x 10. 3 x l°g 10 P + 1.959 X 10~^ x (log-Q p) 2 (2) 

There is good agreement in trends between the curve fit to Sivier' s prediction 
and the present results although the present onset, points fall closer to the 

vapor-pressure, curve than predicted.. If Sivier "s calculations are accurate, the 
results would imply that effects are occurring earlier than predicted for 
homogeneous nucleation. Again, heterogeneous nucieation would be consistent 
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with the data. However, many of the = 0.85 and 0.95 onset points do 
approach Sivier's onset line for homogeneous nucleation. If heterogeneous 
nucleation is indeed occurring then significantly more supercooling may not he 
possible for the M = O .85 and 0.95 tests even if the impurity sites were 
eliminated and homogeneous nucleation took place. 

The present results can also he compared to previous experimental results. 
Goglia in reference 15 has taken data that shows nitrogen onset in the same 
static pressure and temperature range as considered in the present study. 

When Goglia' s data and the present experimental data are plotted together, 
as in figure 37 » the agreement is especially good considering the fact that the 
present data consisted of subsonic and supersonic onset points over a 0.137-m 
airfoil and Goglia' s work consisted of rapid expansions in a 0.1l4-m long nozzle ar 
onset at Mach numbers from 2.0 to 3.25. Although one would normally expect 
more supercooling with the larger velocity gradients in the nozzle, it does 
not show up in the data comparison. In fact, the trend of Goglia* s data falls 
close to the present M to = O'. 75 data but seems to have somewhat less super- 
cooling than was experienced during the present M = 0.85 or O .95 data. Goglia 
states in his conclusion that the condensation he detected 

appeared to be the result of the self-nucleation, or homogeneous condensation, 
of the nitrogen test gas. However, the gas used by Goglia was only 95 percent 
pure nitrogen and contained a 5 percent impurity of oxygen. Even though this 
fact- was well recognized and investigated by Goglia, there may he some effects 
of heterogeneous nucleation on his onset points. (Arthur and Nagamatsu in 
reference 6 detected a mild decrease in supercooling of about 1 to 2 K for a 
5 percent oxygen impurity. ) 
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Whether or not the condensation in the present investigation was the result 
of heterogeneous or homogeneous nucleation, it appears that the magnitude of the 
supercooling is comparable to what others have seen or predicted in the same 
pressure range. 


Benefits Realized 

One of the primary motivations for the present experiments was to determine 
the increase in Reynolds number capability that might be possible when testing 
an airfoil with supercooled flow. This increase was determined for the HACA 
0012-6U airfoil and figure 38 shows the percent increase realized at the tem- 
perature at which condensation effects were first detected. The percent increase 
: was calculated by (l) dividing R corresponding to p and T at onset 
by R corresponding to p, at onset but with T being that which just allows 

C U 

local saturation, (2) subtracting one from the resulting ratio, and (3) multiplying 

by one hundred. As can be seen in figure 38, the M = 0.75 test realized the 

smallest percent gain in capability at 7$» while the = 0.95 test realized 

the largest at 19%. The M ra = 0.75 and 0.85 tests appear to show a decrease in 

gain with increasing p , while a trend for the M =0.95 test cannot 

easily be predicted because of -the lack of a high pressure path data point. 

Error bars shown are calculated from the uncertainty in the value of T for the 

"t 

onset of condensation effects over the airfoil. 

The second motivation for the present experiment was to determine percent 
decreases in drive fan power required and in the amount of liquid nitrogen (LKf2) 
injected when unit Reynolds number is. held constant and the supercooling is 
utilized to decrease both p^ and T^. Since the energy required for producing 
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the amount of M2 injected is much greater than the energy used hy the drive fan, 
reductions in direct operating costs while at test conditions will he much 
closer to the percent decrease in LN2 required than to the decrease in drive 
fan power required. 

The decrease in drive fan power realized is shown in figure 39* The 

trends are, of course, similar to the increase in Reynolds number capability 

shown in figure 38. For M = 01 85 the decrease in drive power required varied 

from 21$ at low pressure to about 19$ at the high pressure end. Decreases of 

23 and 25$ were found for the two low pressure paths in the » 0.95 test. 

The decrease in LN2 required to absorb drive fan heat is shown in figure 1*0. 

For an insulated tunnel such as the Langley 0.3-m tunnel, the injected LN2 

required to absorb the heat flow into .the tunnel through the tunnel shell 
insulation is relatively constant over the temperature range of the present tests 

and is only on the order of 1 percent or less of the amount required to absorb 

the heat of the drive fan. Consequently, heat leakage through the shell has 

no effect on figure 40 and has not' been included. For the =* 0.75 test, the 

reduction in LN2 required varied from 11$ down to 9$ at ‘the higher pressures. 

For' = O. 85 , the reduction varied from 17$ at the lower pressures to 15$ 

at the higher pressures. The M m - 0,95 test had a reduction high of 20$ 

and a low of 18 $. 
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CONCLUDING REMARKS 


An experimental program was undertaken to determine the onset of conden- 
sation effects in the pressure distribution about a 0.137-m NACA 0012-64 
airfoil tested in the Langley 0.3-m transonic cryogenic tunnel at free-stream 
Mach numbers of 0.75 9 0.85, and 0.95 over a total pressure range from Is 2 tq 5*0 
atmospheres. Orifice pressure measurements were recorded at 5$ chord interval's 
and were analyzed to determine the total temperature at which deviations first 
occurred in the orifice pressure coefficient. 

The onset of condensation effects were found to correlate more with the 
amount of supercooling in the free-stream than with the supercooling in the 
region of maximum local Mach number over the airfoil. Effects in the pressure 
distribution over the airfoil were generally seen to appear over its entire 
length at nearly the same total temperature. Both observations suggest the 
possibility of heterogeneous nucleation occurring in the free-stream. 

Comparisons of the present onset results are made to calculations by Sivier 
and data by Goglia. 

The potential operational benefits of the supercooling realized were 
analyzed in terms of increased Reynolds number capability at a given tunnel 
total pressure, reduced drive fan power if Reynolds number is held constant, 
and reduced liquid nitrogen consumption if Reynolds number is again constant. 
Depending on total pressure and free-stream Mach number, these three benefits 
are found to respectively vary from 7 to 19$, 11 to 25$, and 9 to 20$. Since 
operating costs at typical test conditions are approximately 10$ due to fan 
drive power and 90 $ hue to liquid nitrogen injection, reduction in operating 
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costs will be very close to the percentage reduction in liquid nitrogen con- 
sumption. 

Although the present experimental data form a preliminary basis for 
predicting minimum operating temperatures for nitrogen-gas wind tunnels, it is 
just one of several steps that are needed for a more complete understanding' of 
the problem. More work is needed in the following areas: tests of other airfoil 

pressure distributions to determine the sensitivity of the onset of condensation 
effects to different local velocities and gradients, experiments designed to 
determine conclusively if condensation is caused by homogeneous or heterogeneous 
nucleation, studies to determine scaling effects so that the present results may 
be extended to other sizes of tunnels, and studies to extend the results to 
higher total pressures. 
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APPENDIX A 


DATA ANALYSIS PROCEDURE 

Because of the importance of the data reduction technique upon the detection 

of the onset of condensation effects $ this appendix is included in order to 

describe' in detail the correction procedure and to explain the procedure 

for determining airfoil onset T^. Examples of the various steps are included 

for the If = O .85 and R = 16 million path. 

CO c 

The correction procedure was employed in order to take out as much 

Mach number effects as possible when comparing C ' s along a path of constant 

3? 

and R c « While was nominally constant along a path, many pressure 

distributions were recorded at values of M that varied from the nominal value 

00 

by ±.010v To determine if a M ra dependence did exist at a given x/d location j 
a linear regression technique was applied to condensation-free values of 
plotted as a function of M . In this manner, if there was a correlation 
between C' and M , All of the data could be corrected to a single value of 

p CO 

M to . If there was little correlation between C and M^, then a weak dependence 

P 

was assumed to exist and no correction was assumed necessary. Thus all the 
data along a given path of constant M and R could be utilized.- 

OO Q 

The first step of the correction procedure is to apply a linear 

regression technique to values of C at a given x/c as a function of M . 

P . 00 

Using standard definitions as found for example in reference 1 6 , the slope 
of the linear regression fit would be given by 
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The correlation coefficient, r, which is an indication of goodness of fit is 
given by 


r = 


n 

2 ^ 
1=1 i 


Mj(C - C ) 
p i p 


(AU) 


n 

2 (M 

OO 

i=l i 


V 2 .*, <%. - */ 

i=l l 


where if r = 0 there is no correlation and if |r| = 1 then there is perfect 
correlation. For the present analysis if |r| >0.5, then the correlation is 
considered strong enough to use m in correcting a given to another 

value. If [r| < 0.5, the dependence between C and M is considered weak, 

p 00 

and m is assumed to be zero. Three examples of various values of r are 

shown in figure ^1 for the M ra = 0.85 and R c = l6 million path. Figure ^l(a) 

shows the value of for x/c = 0.00 (the orifice at the leading edge of 

the airfoil) as a function of M^. In this case r = 0.15 and, indeed, there 

seems to be little correlation between M and C . Consequently, no attempt 

• P 
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is made to manipulate or correct this data for small differences in M^. 

Figure 4l(b) shows the values of C for x/c = 0.40 as a function of M^, 

In this case r = 0.49, and so this linear curve fit barely fails to meet the 

correlation test. Figure 4l(c) illustrates the values of C for x/c = 0.80 

that lead to a value of r equal to 0.92. As seen in figure 4l(c), there is 

good correlation between C and M , and differences in C due to differences 

P P 

in M can be removed by correcting values of C to a given M by setting 


C = C - m • (M - M ) 

P P CO CO 


where C is the corrected value of C and m is the slope of the linear 
P 0 P 

fit. Again, m is assumed to be equal to zero if |r| < 0.5. 

Onc'e the linear regression technique has been applied to the values of 

■for each orifice along the path, a set of intermediate graphs are drawn before 

the airfoil onset T, is determined. For each orifice a value of AC is 

t p 

plotted as a function of where 


AC = C - C 


P P 

■^o ^ave 


In equation (A6), C is corrected with respect to M and C 


is an 


average of condensation -free values of C for that orifice. (Condensation- 

1 P 
•^o 

free values were determined by establishing crude onset values of without 

the benefit of correction procedures and then adding at least a 2 K temperature 


buffer. Final onset values of T, were then checked to insure that C 

t P 

ave 

was indeed free from condensation effects.) An example of a AC vs. T, 

p t 
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with and without correction is shown in figure 42 for x/c = 0.80 and where 
r = 0.92. Figure 42(a) shows the uncorrected data while figure 42(b) shows the 
values of C after being corrected. As can be seen, considerable "order" is 

ir 

added to the data by the correction process and this "order" permits easier 
identification of the T^_ at which C at this x/c locatdoft begins to 
depart from its unaffected average value. Hot always, however, was |r| > .5 
and so this simplification was not always possible. Graphs similar to fig- 
ures 42(a) or 42(b) are then drawn for each orifice location over the airfoil. 

These graphs are used to determine at what T. each of the orifices first 

t 

experience systematic departure from its unaffected average value of C 

^0 

These orifice onset temperatures are then plotted as a function of x/c 

as was shown in figure 20. These orifice onset T vs. x/c plots are the 

key analysis graphs because they provide an overview of the effects on the 

airfoil as a whole. There is experimental uncertainty in the onset point for 

each orifice and consequently it is helpful to see all of the orifice onset 

T-f. ’ s together so that a particularly low or high value can be checked for 

possible error. The airfoil onset T^_ is not taken to be an average of 

orifice onset T 1 s, but rather is taken to be the warmest orifice onset T 
■ « t 

that appears systematic. Uncertainty bars are normally drawn to encompass 
the highest orifice onset T^, whether the deviation appears systematic or not. 

The data analysis procedure described herein was more sensitive to 
condensation onset than previous analyses described in references 17,- 18 and 
19* In references 17 and l8 pressure distributions were superimposed to detect 
differences - no correction was made to C for differences in M and. 

p 00 5 

consequently, much of the data could not be used. Reference 19 utilized the 
AC? vs. Tj. format at 0, 25', 50, and 75 percent chord locations but still did 
not include the Mach correction technique. 
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APPENDIX B 


ERROR ANALYSIS FOR AC VS. T -GRAPHS 

P t 


During the present experiment, fluctuations in wind tunnel conditions were 

the prime source of error during the 50 second period required to record the data. 
Tunnel conditions were observed to fluctuate with standard 


deviations of the following approximate magnitudes: 0 of + .005 atm (varied with M ot 

t 


R c h of ±.003; and 0^ of ±0.5 K. The manufacturer’s stated uncertainty in 

00 -fc 

the static pressure measurements, Ap, was 0.5# of full scale. A representative 

standard deviation for AC due to these fluctuations is calculated herein 

P 

using normal error propagation techniques, see reference 20 for example. 

As defined in equation (A6), 



•ave 


(A6) 


in which C is the value of the pressure coefficient, corrected to some 
p o 

value of free-stream Mach number, M^, from data taken at a free-stream 

© 

Mach number of M ot . The correction is defined as 



(A5) 


where m is the slope of the linear regression fit as defined in equation (Al). 

If sufficient correlation, r, as defined in equation (aU), is not found, m 

is assumed to be equal zero. Since C is just the average of the corrected 

P ave 

condensation-free data, the variance of C can be approximated by 

P ave 
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;2 


w 


(Bl) 


aye 


where - If would be the number of, samples used in the average. The variance of 

C can be written from equation (A5) as 
-^o 


2 2 2 
°C = ’ a C +a F 


(B2) 


where F is defined as 


F = m(M - M ) 


OO CO 


(B3) 


Assuming there is much more relative error in (M - ) than m, the variance 

CO CO ' 


in F will be primarily due to the variance in the difference in Mach numbers. 

Since is an arbitrarily defined constant, we may write' 

o . 


2 _ 2 2 
F ~ n 


(b4) 


Consequently, equations (A5 ) s (a 6), (Bl), (B2), and (B3) lead to 


_2 + l w _2 . 2 _2 ^ 

°AC = (_ T — c + ™% ) 
p p 00 


(B5) 


It now remains to discuss a c , the variance associated with error in the 

■ P 

values of C before Mach correction. The value of C is defined as 
P P 


P - P 0 

C = 

P %> 


OF POOR QUALII* 


(B6) 
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Though Cp is strictly a function of B and M as Schlichting explains in 

reference 21, the calculated value of C from, an experiment may he sensitive 

to errors or fluctuations in p» P^.> T^_, and M^. Tn other words, one must 

distinguish between actual C , C (act), and calculated C , C (cal). During 

P P P P 

the present experiment, the prime difference between (act) and (cal) 

was lack of simultaneity - tunnel conditions were recorded at some time t and 

the static value of p at an orifice may have been recorded at a time up to 50 

seconds later. Calculation of the impact of all the fluctuations upon the 

possible error in C (cal) will be presented below except for the dependence 

P 

of C (cal) on T^_. Since only enters the dependence of (cal) through 

R, and since the R dependence was weak, no error in due to was 

calculated. An equation for the variance of can be written as 


ac 


3C 


3C 


(cal) 

P 


(RT^) cl + (g^) 0 ,, + (t^t) a «' 


'3p + 


M_ 


' 9Ap Ap 


;(B7) 


The next three sub-sections describe the calculation of the partial derivatives. 


Solving for 


3C p (cal) 

9 P + 


'V R > 7 


Equation (b 6) can be rewritten as 


C (act) = 
P 



P 


t 


(B8) 
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Since M and R are held constant during the process, then C (act), — 


P-t 

yields . 


and — all remain constant. Thus for a p fluctuation, equation (b 8) 
p t t 


R = 


= Constant 


or that 

ie_i x = 

8pjM , R, - 

°° c -^t 

Now, turning hack to C^(cal ) and remembering that , Poo and thus 
are all recorded at some time t^ and then p is recorded later at some 
equation (B8) can be written as 


so that 



9C (cal) 

_ E 1 x 

3p, 'M . R, - 

t 00 c 


1 a P | 

\ 9p t M co* R ’ C 


or with equation (BIO) and dropping the subscripts "1" 


9C (cal) 

_E 1 x = J£— 

3p t ! M ot , R, - 


(B9) 


(BIO) 


t. 


(Bll) 


(B12) 


(B13) 
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Assuming the isentropic relation for a perfect gas with its ratio of specific 
heats equal to 1.% equation (B13) can he written as 


ac.(cal) 

P 

3 p t~ 




+ 1 ) 


(BlH) 


r. X 

This partial derivative is found when the linear regression step in the 
Mach correction technique is performed as outlined in appendix A. The 
value of 3C (cal) is equal to the value of the slope, m. 


Solving for 3 ^p^ cal ) 

§m“ 


30 (cal) 

p = m 

3M m 

00 

because m is just the linear regression fit to calculated values of C 


(B15') 


as a function of M ot . As before, when the correlation coefficient, r, 
does not reflect good correlation, then 3C^(cal), like m, will be assumed 

Wr 

00 

equal to zero. 


3C t 

Solving for p j 
3Ap 1 

■ This is just the simple case where C (cal) is in error because of 
instrument error in p. That is, 

p(cal) = p(act) + Ap 


(Bl6) 


so 

^2 I x “ 1 

3Ap |P t ,M oj , R,- 
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Consequently, 

3C p (cal) = 1 

9 A P <Jco 

or in terms of M m and p^ for the isentropic flow of a perfect gas with a 
ratio of specific heats equal to 1.4 


(B17) 


SC (cal) . (1 + *) 

p 5 


^, 3 - 5 


SAp 




(B18) 


At this point, all of the partial derivatives in equation 

(B7) are known and can "be combined with the observed or known values of 

a , a„ , and a. The value of 0, was first taken to be the 
P t V Ap. A P 

manufacturer's value of 0.5 percent of full scale, but this value gave a 

contribution to the error bar in C that was on the order of 10 times 

P 

the actual fluctuations in C . Apparently the two pressure transducers 

p 

performed- with much better precision than expected. In fact, instead 
of the expected respective uncertainties of 0.0170 and 0.0085 atm for 
the M = 0.75, 0.95 and the = 0.85 tests, the uncertainties were on the 
order of .0004 and .0002 atm, respectively. These low magnitudes in % 

made the contribution of the instrument error to the total error negligible. 

Assuming isentropic relations for a perfect gas with the ratio of 
specific heats equal to 1.4 and assuming <J^ is negligible, equations 
(35), (B7), (Bl4), and (B15) can be combined to give 


2 2 


= ( * . + . I) a p + 2 m cr } 

n W p .7*C p 


(B20) 


The error bars used in the AC vs. T, graphs in the data analysis procedure 

P t 

AC * 


were calculated with the above value of a 
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Figure 1.- Effect of temperature reduction for = 1.0, constant P. 
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and constant reference length. 
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c = 0.137 m 


Figure 3.- Two dimensional NACA 0012-64 airfoil has orifices 
spaced at 5 % chord locations. Trailing edge orifice 
was added for tests with M =0.75 and 0.95. 
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Figure 7.- Paths of constant r q for the.M^ = 0.95 test. 
(R c a 40 X 10 6 path taken with M m jb 0.9,3.) 







•Figure 11.- Orifice onset T t * s and resulting airfoil onset T t 
for the R c = 15 x 10 6 path of the M w = 0.75 test. 
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Figure 12.- Orifice onset T t »s and resulting airfoil onset T 
. • f° r R c = 24 x 10® path of the = 0.75 test/ 
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Figure 14.- Onset of condensation effects for M =0.75 test. 
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Figure 16.- Average unaffected pressure distribution for the 

R = 26 X 10 6 path of the M =0.85 test, 
c x °° 






Figure 18.- Average unaffected 1 pressure distribution for the 
R c ” 38 X 10 path of the = 0.85 test. 
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Figure 19.- Average unaffected pressure distribution fpr the 
R = 42 X !0 6 path of the M - 0.85 test. 
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Figure 20.- 


Orifice onset T t ’ s and resulting airfoil onset 
for. the R = 16 x 10 6 patfi of the M =0.85 test. 
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Figure 21.- Orifice onset T t ’ s and resulting airfoil onset T t 
for the R c = 26 x 10 6 of the M ro = 0.85 test. 
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Figure 22. - Orifice onset ' s and resulting airfoil onset 
R c = 34 X 10 6 path of the M w = 0.85 test. 
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Figure 23.- Orifice onset T t T s and resulting airfoil onset T t for the 
R = 38 X 10® path of the M =0.85 test. 
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Figure 24.- Orifice onset 
R = 42 X 10 6 
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* s and resulting airfoil onset for the 
path of the M =0.85 test. 
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Figure 26.- Average unaffected pressure distribut 
R = 16 x 10 6 path of the M =0.95 
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Figure 29.- Orifice onset 
R = 16 X 10 6 


T t f s and resulting airfoil onset for the 
path of the = 0.95 test. 
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Figure 30.- Orifice onset T t ’ s and resulting airfoil onset T. for the 
R c = 27 X 10 6 path of the = 0.95 test. * 
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R c =• 40 x 10 path of the !M = 0.95 test. 




Figure 32.— Ac vs, T for M — 0,93 and R = 39*8 million 

u w C 

Unaffected average value of C ■ =' -0.035. 
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Figure 33.- Onset of condensation effects for the M = 0.95 test. 
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Figure 34.- Supercooling for all tests. 
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(b) Static values based on M = M . 

Figure 35.- Concluded. 
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Goglia 1 s data from reference 18. 






Figure 39.- Percent decrease in drive fan power at onset. 
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•Figure 40.- Percent decrease in LHj required to absorb drive 

fan heat. - 






(a) Before Mach correction. 

Figure k2.~ Improvement in graphs of ACp ts. T with Mach 

correction for M = 0 . 85 and R = l6 million'. 
00 . C 

Unaffected average value of Cp =? -0.0li2. 
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